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Abstract: We analyze a continuous spin Gaussian model on a toroidal triangular lattice 
with periods Lq and L\ where the spins carry a representation of the fundamental group of 
the torus labeled by phases uq and u\. We find the exact finite size and lattice corrections, 
to the partition function Z, for arbitrary mass m and phases U{. Summing Z -1 / 2 over 
phases gives the corresponding result for the Ising model. The limits m —>■ and U{ — > 
do not commute. With m = the model exhibits a vortex critical phase when at least 
one of the U{ is non-zero. In the continuum or scaling limit, for arbitrary m, the finite 
size corrections to — In Z are modular invariant and for the critical phase are given by 
elliptic theta functions. In the cylinder limit L\ —>■ oo the "cylinder charge" c(uq, vr?L§) 
is a non-monotonic function of m that ranges from 2(1 + 6uq(uq — 1)) for m = to zero 
for m — > oo. 

PACS numbers: 05.40.-fj, 05.50.+q, 64.60.-i, 11.15.Ha, 05.70.Fh, 75.10.Hk 



Finite size effects are an intrinsic feature of laboratory experiments which probe the neigh- 
bourhood of a continuous phase transition. The scaling properties of these corrections to 
the infinite (bulk) system behaviour play an increasingly important role in computer sim- 
ulations and our theoretical understanding of the critical regime of statistical systems [1] . 
The interplay of these effects with other aspects of the system can give rise to crossover 
from one characteristic behaviour to another [2] . It is possible to make real progress in the 
detailed analysis of such situations in two dimensions and for this reason two dimensional 
models have attracted much interest in recent years. Many of their properties are express- 
ible in terms of those of generalized Gaussian models [3] and with minor modification the 
results of this note can be adapted to a wide class of other models. 

Pure finite size effects can be isolated from those due to the presence of a boundary by 
ensuring that the connectivity of the underlying lattice is such that it has no boundary. 
One can avoid other complications, such as local curvature of the lattice, while retaining 
the finite size effects, by considering a flat torus. Furthermore, in two dimensions, the only 
zero curvature finite volume manifold without boundary is the flat torus. This torus can 
be conveniently thought of as a parallelogram with opposite sides identified. We take the 
sides to be of length Lq and L\, with Lq at an angle 9 to L\, cf. fig. 1. Simple geometry 
then means that the point (x, y) is thereby identified with the point (x + aLi cos 9 + bLQ, y + 
aL\ sin 9) where a and b are integers. 

We consider a triangular lattice composed of similar triangles, pairs of which form 
parallelograms, cf. fig. 1. The basic triangles have two sides of lengths oq and a\ with an 
angle 9 between them. The complete lattice forms our torus, T 2 , and consists of KqK\ 
sites and 2KqK\ triangles, forming a parallelogram of sides Lq = Kqoq and L\ = K\a\. 

We take a complex continuous spin variable ip(x, y) on T 2 but will not demand that ip be 
periodic, though the underlying lattice is, rather it acquires a phase on being transported 
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around one of the cycles of the torus. In general we require that 

ip(x + cos 9 + bL Q ,y + aL 1 sin 9) = e 2m ( aUl+bu °^(x,y), a, b G Z (1) 



For uq and u\ neither integer or half integer the spin variable is necessarily either complex or 
real with two components: From the two component viewpoint these boundary conditions 
can be thought of as corresponding to vortices winding around the periods of the torus. 
From a mathematical standpoint ip is a section of a bundle C over the torus T 2 . 

Labeling the lattice sites by (/co, k\) = k = Icq + {k\ — 1)Kq with k$ = 1, . . . , the 
energy of a configuration is given by 

S L [T\ <p\<p\ = \Y, v^*(k) (A(k, k') + m 2 d Kk ,) V (k') (2) 
kk' 

where ^fg = aQa\ sin 9 and S^ ^i is the Kronecker delta. A is a KqK\ x KqKi symmetric 
matrix, all entries of which are determined by nearest neighbour interactions, so that each 
spin interacts with six neighbouring spins. Explicitly the only non-zero elements are those 
given below (and their transposes): we have 

h), (k + 1, ki)} = -a = — ^y- I \ - ) 

sin z 9 \ag agai / 

A{(k , h), (k Q , k 1 + l)} = -(3 = - * (I - — ) 

sm z 9 \ af agai / 

(3) 

A{(k + 1, fci), (fco, h + 1)} = - 7 = - * (—) 

sm z 9 \ a a l / 

A{(/co 5 k\), (/cq, ^i)} = 2a, where a = a + (3 + 7 

We actually solve the model for general nearest neighbour couplings, a, (3 and 7, in which 
case g = (aj3 + (3^y + 7a) -1 . Despite the simplicity of the model we will see that it has a 
surprisingly rich structure. 
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The eigenfunctions of A that satisfy (1) are the discrete Fourier basis 
1 



(k ,ki) 



exp 



2ttz { (m + u o)j^ + (n + u l)j^ 



with 



orthogonality of which is given by 



K K x 

e mn(^0? ^l) e m'n'(^0) ^l) = ^m,m'^n,7 

k =l fc =l 



/ m = 1, . . .K 
\n = l,...K 1 

(4) 



(5) 



while the corresponding sum over m and n gives completeness. The eigenvalues X no ni of 
the matrix |(A + m 2 ) are then 

A no m = 5-acos(x no ) -/3cos(x ni ) -7cos(x no - x ni ) 



2n(rij + ui) . n ^ . r m 

where; .r„ ( = — , / (). 1 and o a ■ — — 



(6) 



Ki 



We are interested in the partition function 
Z(T 2 ,£,m) = J [Jjd^*d^ 
The free energy is then given by F = k B TW where W = — In Z 



e k b 



{K ,K{) , 
(no.m) Lv^nom 



(7) 



To perform the sums we note the rearrangement of the eigenvalues 

An n! = 5 ol cos a; no - |/? no | cos(a; ni - 6 no ) where (3 no = (3 + -fe tXn o = \p no \e l6n ° (8) 
and use the basic non-trivial identity [4] 



K-l 
n=0 



COSXr 



10 7T 



2 2 
Z — COS V 



+ In 



1- (z- y/z 2 -l^ K e 27riu 



(9) 



where z > 1 and real. In any case the sum over n\ followed by that over uq yields: 



W = K K 1 W B + W F 



(10) 



where 
W B = 



71 dui r du 2 

-7T 2yr J-ir 2yr 



ln 



nk B T 



{5 — acos(z/i) — ,9cos(z/2) — 7cos(z/i — z^)} 



(11) 



and 



Wp =K X 

e=± 



K 

1 - [Qe ~ \lqt - 1 ) e 



(12) 



here [x] denotes the integer part of x. The definitions of the various variables in (12) above 
are as follows: 



q±=q±\ q z -p q 



aS + /?7 cos 2 v 5 2 - (/3 2 + 7 2 ) cos 2 v 



or 



P 



CT 



(13) 



and 



Vn 



-ln 



^no - \J z n ~ 



+ i0. 



n 1 



where 



Zn 



5 — a cos Xn 



2|/3n | 

Note that Wb gives the free energy per lattice site in the thermodynamic limit, i.e 

W 



(14) 



lim 



W B 



(15) 



hence, due to the fact that we are dealing with a finite rather than an infinite lattice, Wp 
gives the complete finite size correction to the bulk lattice behaviour. 

If instead we perform the sum over no before that over n\ we find the alternative 
composition W = KqK\Wb + Wp. The two apparently different expressions Wp and Wp 
are, obviously, equal; they can be transformed into one another by the interchange of a 
with (3 and the subscripts and 1. 

The limit of interest to us is the continuum or scaling limit. This is a constrained 
thermodynamic limit achieved by taking Kq, K\ — > oo while keeping fixed Li = K^ai, 9, 
m 2 and the ratio k = The asymptotic form of Wb in this scaling limit is given by 



T/to 2 T/to 2 to 2 !/ 

WbKqK! = KqK^b - ^— {HKoKl] ~ M - ^(M^l " 1) + 



4n 



47T 



(16) 



Here is the value of Wq at m = and 



PIT 

P = p(a, (3,y) = / 
JO 



dv 



1 



1 



sin v 



\J\ + (7a 2 sin 2 f ? ' 



iln[vW + 7)] (17) 



Despite its appearance, p(a, j3, 7) is symmetric under interchange of a and (3. Both 
and p depend on Lq, Li, 9 and k, or equivalently they depend on the geometry of the 
lattice triangle. 

We find the limiting behaviour of Wp is given by(«) 



00 



tit} , m V. \ - 
r F = -^c{u , )+ V In 

T 1 f * 
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n=— 00 



^ -27TT 1 y / '(n+uo) 2 +^i^+27r J {u 1 -To(n+u )} 



where 



fc 



^0 = 



n = 



fey 



vW + 7l' " /^ + 7 
The function c(tt, x) that appears in (18) is given by 



and V = K§K\^fg 



00 



dp 
2k 



In 



2 _ e -y/p 2 +x+2iriu 



c(u, x)n 
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(18) 



(19) 



(20) 



Corrections to (18) which vanish in the limit are also easily obtainable from (12). In 
terms of the torus geometry described above we have tq = cos 9, t\ = sin 9 and 
V = LqLi sin6>. 

Certain special cases of (18) arise in the case of an Ising model on a torus (see [5] 
where the case tq = was studied). This can be understood from the equivalence of the 
Ising model to a dimer model on a decorated lattice where similar determinants arise [6]. 
From this equivalence, exhibiting the phase dependence of Wp by Wp(uQ, u\), with a, (3, 
7 and S obtained from [7], we have 

Z Iaing = l^-W 1 ™ 9 | Te lT^ F (0,0) + e \W F {Q,\) + e \W F {\,Q) + e |^(I,I) j (21) 



( a ) In extracting the limit, the expansions lng_ = — -^(/3 + 7) 2 (^ 2 + ^pj;) + ••• and 



1 



9(P + 7)1/3 + 71 



4. i — — — x„ n + ■ ■ ■ for small v and m 2 are useful. 



1/3 + 71 /3 + 7 



for ferromagnetic couplings, with + referring to T < T c and — to T > T c . In the scaling 
limit Wp — > Tp. Our results therefore incorporate the complete lattice and finite size 
corrections for the Ising model on a triangular lattice. With a similar equivalence our 
results can easily be translated to give the general result for other models. 

Tp is invariant under the transformations 

(i) UQ I— > UQ U\ I — > UQ -\- U\ TQ \ — > 7~0 + 1 T\ I — ► T\ 

TQ T\ 



(ii) UQ I — ► U\ U\ I — > —UQ TQ 



|2 1 I 1 9 

rr rr 



where t = tq + %t\. Invariance under (i) is obvious from (18) and records the action of the 
generator 'T' of SL(2, Z)/Z2 while (ii) interchanges Tp with Tp this latter being obtained 
by performing the sums leading to Tp in the opposite order. These two expressions are 
equal and (ii) records the action of the generator '5" of the modular group SL(2, Z)/Z2- 
The volume V is invariant under the action these generators. But S and T generate the 
whole modular group SL(2, Z)/Z2 and so Tp and V are invariant under the action of any 
element of SL(2, Z)/Z2- Hence we can conclude that in the scaling limit the complete finite 
size corrections to the bulk free energy are invariant under the entire modular group. 

If M denotes an arbitrary element of the modular group, under which r undergoes the 
well known transformation r i— > (ar + b)/ (cr + d), we have that Tp is invariant under the 
replacement 

u\ t— > ui[M] = au\ + buQ uq h- > uo[M] = cu\ + duQ 

(22) 

Li ^ Li[M] = aL x + bL Q e' td Lq h-> L [M] = cLie** + rfL 

where 

M=(^ c jjeSL(2,Z)/Z 2 , a,b,c,deZ, ab-cd = l (23) 

We note that the sum over the four terms in (21) gives the finite size correction Z^ ng 
as a modular invariant function of tq, t\ and m?V . It is clear from our construction that 
summing over the phases U{ which form an orbit of SL{2, Z)/Z2 on the space of phases 

7 



(the Picard Variety [8]) will allow one to construct similar phase independent modular 
invariant partition functions, such as arise in other conformal field theories cf. [9] and 
references therein. 

Taking m — > in (18) we find the critical phase limit of Tp is given by 



I> = - —c(u ) + ln 1 ~ e 27r ^ T l n l- e «K-™o)} 



oo 



In 



6 

n=— oo 



n + UQ 
\n + UQ\ 



(24) 



7](t) 

where c(u) = 2{1 — 6u(l — u)}, t?i is a Jacobi theta function and r\ is the Dedekind eta 
function — the modular transformation properties of d\ and r\ are well known and allow an 
independent check of the modular invariance in this limit. 

A further limit of interest is the geometric one obtained by taking L\ — > oo; it corre- 
sponds to a cylindrical geometry. In this limit Yp/V reduces to the expression ^cylinder 
where 

7 cylmder = --^c(u ,m 2 L 2 ) (25) 

(The analogous limit, where Lq — > oo, replaces uq and Lq by u\ and L\ respectively.) 
When uq and m are both zero the "cylinder charge" c(uq, m 2 Lg) reduces to the central 
charge of the model; in our case c = 2. The cylinder charge should not be confused 
with the Zamolodchikov c-function [10] for this model, the two functions have different 
dependencies on m, the latter for example being a monotonic function of m. We plot 
c{u,x) in fig. 2 for different values of u as a function of x. It exhibits the crossover to 
c = at large values of x, and is clearly not a monotonic function of x. From (21) we see 
that 

(26) 

Comparison with (25) means that the cylinder charge for the Ising model is —^c(^,x), 
which for x = gives the usual central charge c = \ [11]. 
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The limits m — > and uq, u\ — > for r_p do not commute. To see this we expand IV 
(18) for small m and uq, u\ obtaining 



and this expression clearly can tend to the distinct logarithmically singular expressions 
\n\u\ — rtto| 2 + 21n|?7(r)| 2 and ln[rim 2 y] + 2 In |?7(t) | 2 depending on the order in which 
the limits are taken. However both limits and indeed (27) itself are modular invariant. 

The continuum version of the model studied here is described by the Hamiltonian 



It is a simple task to obtain the corresponding continuum partition function by the method 



an arbitrary undetermined scale. Analogues of this continuum model were discussed in 
[9] where the phases were taken to be rational numbers. Their expressions unfortunately 
contain errors and they did not succeed in separating the bulk and finite size contributions. 

The model we have studied above serves as a useful starting point for a perturbative 
treatment of the approach to the Kosterlitz-Thouless phase [12] in a |^| 4 theory and thus 
captures the neighbourhood of the critical point of the XY model as the critical point is 
approached from the disordered phase. That the boundary conditions capture the essential 
features of a vortex phase can be seen by considering (1) in the cylinder limit and mapping 
the cylinder to the plane using the conformal map y + = Inz. The charge of the 
vortex at the origin is then given by uq and its existence gives rise to an Aharanov-Bohm 
effect for transport around the origin. 

In summary, the finite size corrections to the free energy are modular invariant. This 
conclusion extends to the entire scaling neighborhood of the critical phase. We use our 




(27) 




(28) 



of ^-function regularisation. This procedure yields the finite size contribution of (18) 
but the bulk term differs significantly being now — ^-(ln[(m/27r / u) 2 ] — 1) where \i is 
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results to give expressions for the complete lattice and finite size corrections for the two 
dimensional Ising model on a triangular lattice via its equivalence to a sum over Pfaffians. 
Modular invariance also extends to models in more than two dimensions when the geometry 
giving rise to finite size effects contains a flat torus. For a three dimensional cylindrical 
geometry with toroidal cross-section the result can be obtained from by replacing m 
with m 2 +(/ 2 and integrating the resulting expression over q. One can understand the origin 
of modular invariance in general as the residual freedom to reparametrize coordinates, in 
the continuum limit, while retaining flat toroidal geometry. 

In the two dimensional case the limiting finite size corrections at the critical phase are 
expressible in terms of classical elliptic functions. Infinitesimally small values of the phases 
Ui lead to logarithmically divergent contributions to the free energy. This implies that the 
free energy needed to create a vortex becomes infinite for an infinitely large lattice. In 
general the model has a surprisingly rich structure of non-commuting limits. For example 
the limits of approaching the critical phase and that of sending the to zero do not 
commute. 

Acknowledgment: We are grateful to Paul Upton for helpful conversations and his 
careful reading of the manuscript. 
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The triangulated torus 




Fig. 1. 



The 'cylinder charge' function c(u, x) for various u 



c (u, x) 




Fig. 2. 



11 



References: 

[1] M. N. Barber, in Phase Transitions and Critical Phenomena, vol.8, eds. C. Domb and 
J. L. Lebowitz (Academic Press, London 1983). Finite Size Scaling CPSC vol.2, ed. 
J. L. Cardy (North Holland, 1988). 

[2] D. O'Connor and C.R. Stephens, Int. J. Mod. Phys. A9, 2805 (1994). 

[3] VI. S. Dotsenko and V. A. Fateev, Nucl. Phys. B240 312 (1984), and B251, 691 
(1985). 

[4] This is a generalization of an identity in A. E. Patrick, J. Stat. Phys. 72, 665 (1993). 
[5] A. E. Ferdinand and M. E. Fisher, Phys. Rev. 185 832 (1969). 

[6] P.W. Kastelyn, J. Math. Phys. 4, 287 (1963). M. E. Fisher, J. Math. Phys. 7, 1776 
(1966). 

[7] J. Stephenson, J. Math. Phys. 8, 1009 (1964). 

[10] A. B. Zamolodchikov, Pis'ma Zh. Eksp. Teor. Fiz. 43, 565 (1986) [JETP Lett. 43, 
730 (1986). 

[11] H. W. J. Blote, J. L. Cardy and M. P. Nightingale, Phys. Rev. Lett. 56, 742 (1986). 
I. Affleck, Phys. Rev. Lett. 56, 746 (1986). 

[8] D. B. Ray and I. M. Singer, Ann. Math. 98, 154 (1973). 

[9] C. Itzykson and J. Drouffe, Statistical Field Theory vol 2, C.U.P. (1989). 

[12] J. M. Kosterlitz and D. J. Thouless, J. Phys. C 6, 1181 (1973). 



12 



